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Abstract The different methods of determination of black holes (BHs) 
masses are presented for three classes of BHs observed in the Universe: 
stellar mass BHs, intermediate mass BHs (IMBHs) and supermassive BHs 
(SBHs). The results of these determinations are briefly reviewed: stellar 
mass BHs are found in the range of about 3 to about 20 Mq , IMBHs in the 
range of a few hundreds to a few tens of thousands Mq (the determinations 
are much less precise for these objects) and SBHs in the range of about 
3 X IO^Mq to about 6 x 10^° Mq. 

Key words: — black holes: mass determination — stars: black holes — 
intermediate mass black holes — supermassive black holes 

1 INTRODUCTION 

Black holes (BHs) observed in the Universe can be classified into three groups: stellar 
mass BHs, intermediate mass BHs (IMBHs) and supermassive BHs (SBHs). The case 
for the existence of IMBHs is still substantially weaker than the case for the two other 
groups (the existence of which is, practically, beyond any dispute). Therefore, the mass 
estimates for these two groups are substantially more precise than for IMBHs. This paper 
is devoted to a brief review of the different techniques used to estimate the masses of 
BHs and to the brief presentation of the results. 

2 STELLAR MASS BLACK HOLES 

As of today, the stellar mass black holes, essentially, are observed only in the X-ray 
binaries (XRBs). Only a few rough estimates are available now for single black hole 
candidates from microlensing events (see the next subsection). 
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2.1 Black Holes in X-Ray Binaries 

At present, 58 black hole candidates are known among compact components of XRBs. 
The masses are determined for 24 of them (see Table 1). Some of these determinations 
belong to the most precise mass estimates ever derived for any black hole. Below, I briefly 
summarize the technics used to obtain these estimates. 

• The mass function 

This is the most important observational parameter used to constrain the mass of the 
compact component. The mass function /(M^) is calculated from the radial velocities of 
the optical companion: 

/(M,) = 1.0385 X W-^K^j,tP Mq (1) 

where Kopt is the scmiamplitudc of the radial velocities of the absorption lines of the 
optical component (in km/scc) and P is the orbital period (in days). The mass function 
is related to the masses of both components by: 

/(M,) = MlsinH/{Mopt + M^f (2) 

where and Mopt are the masses of the compact and the optical components and i 
is the inclination of the orbit. The value of f{Mx) gives an absolute lower limit to the 
compact component mass. To obtain a more precise value (not just the lower limit), we 
have to estimate the mass ratio q = Mopt/M^ and the inclination of the orbit i. 

• The rotational broadening of the absorption lines of the optical compo- 
nent 

From the measurements of the lines, we determine the projected rotational velocity 
at the equator of the optical component Vrotsini. Assuming the corotation of the optical 
component with the orbital motion (for a Roche lobe flUing component, it is a very good 
assumption), we have: 

Vrotsini = 0A6Kopt{q/{l + q))^^^ (3) 

With the help of this equation, one can determine the mass ratio q. The intrinsic width of 
the absorption lines is very small ( ^ few km/sec), while the typical rotational broadening 
is of the order of few tens km/sec and can be. relatively easily, measured. This method 
has been applied, with a substantial success, to many BH systems (sec Orosz, 2003). 

• The radial velocities of the emission lines of the accretion disc 

The shifts of these lines reflect the orbital motion of the compact component and so 
permit us to determine the mass ratio directly: 



Q — K em/ Kopt 



(4) 
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where Kf-m is the semiamplitude of the radial velocities of the emission lines. This 
method of determining the mass ratio is completely independent of the previous one. 
Unfortunately, it is rather uneasy to implement, since the emission lines are very broad 
(~ 2000 km/sec), while the orbital shifts are of the order of few tens km/sec. 

• The amplitude of the eUipsoidal Hght variations 

Due to filling of the Roche lobe, the optical component is tidally distorted and due 
to rotation it exhibits the ellipsoidal light variations (doiiblo sinusoid per orbital period 
if the rotation is synchronous). The amplitudes of these variations in V and I are given 
by: 

= 0.26smV(l + 9) (5) 

M = Q.2AsinH/{l + q) (6) 

As may be seen, the dependence on the mass ratio q is rather weak (for most of the 
BH XRBs q falls in the range 0.05 0.2) and, therefore, the ellipsoidal light variations 
provide us with a valuable information about the inclination i (even, if we do not know 
q). In practice, the procedure is not that simple, since the optical light is usually (even 
in the quiescence) contaminated by the residual contribution from the accretion disc. To 
make the situation worse, this residual contribution is frequently variable. Extraction of 
the true value of the inclination requires, often, a very careful modeling (Froning and 
Robinson 2001, Gclino ct al. 2001). 

• The mass-spectral type relation for the optical component 

The typical optical components of BH XRBs are the lower main sequence stars, 
which satisfy reasonably well the mass-spectral type relation. This relation may be used 
to estimate the mass of the compact component, if the mass ratio is unknown. Since 
the mass of the optical component is usually quite small (below 1 Mq), even substantial 
uncertainty (say factor of two) does not influence dramatically the mass estimate for the 
compact component. 

• High frequency QPOs and X-ray spectra 

Using this method, one can estimate simultaneously the spins and the masses of 
black holes. The spin can be estimated from careful analysis of either continuum X-ray 
spectrum (Shafee et al., 2006; Davis et al., 2006; McClintock et al., 2006) or spectral X- 
ray Hues (Miller et al., 2002; Miller et al., 2004; Miller et al., 2005; Miniutti et al., 2004). 
The first approach gives, at present, the more reliable results. Once the spin is known, 
we can use the observed high frequency QPOs to derive the mass of a given black hole. 
To do so, one has to apply one of the theories of high frequency QPOs in the systems 
containing accreting black holes. The most successful one seems to be, at the moment, 
the parametric epicyclic resonance theory (Abramowicz & Kluzniak 2001, Abramowicz 
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et al. 2004, Kluzniak et al. 2004, Lee et al. 2004, Torok ct al. 2005). For a brief summary 
of this theory, the reader is referred to e.g. Ziolkowski (2007). 

The results of mass estimates carried out with the help of all the technics described 
above are given in Table 1. 



Tab. 1 — Masses of Black Holes in X-Ray Binaries 



Name 


Porb 


Opt. Sp 


X-R 


C 


Mbh/ Mq 


Rcf 


Cyg X - 1 




09.7 lab 


pers 




20 ± 5 


1 


LMC X - 3 


l''70 


B3 V 


pers 




6-^9 




LMC X - 1 


4''22 


07-9 III 


pers 




4 10 




SS 433 


13"^ 1 


~ A7 lb 


pers 


mQ 


4.4± 0.8 


2 


LS 5039 


3'^906 


07f V 


pers 


mQ 


2.7-^ 5.0 


3 


XTE J1819 - 254 


2''817 


B9 III 


T 


mQ 


6.8 ^ 7.4 




GX 339 - 4 


1''756 


F8-G2 III 


RT 


mQ 


6 


4 


GRO J0422 + 32 


5''09 


M2 V 


T 




4 ± 1 




A 0620 - 00 


7'' 75 


K4 V 


RT 




11 ± 2 




2S 0921 - 630 


9^*01 


KO III 


pers 




1.7 4.3 


5,6 


GRS 1009 - 45 


6'' 96 


K8 V 


T 




4.4 ^ 4.7 




XTE J1118 + 480 


4'n 


K7-M0 V 


T 




8.5± 0.6 


7 


GS 1124 - 684 


10''4 


KO-5 V 


T 




7.0 ± 0.6 




GS 1354 - 645 


2^*54 


GO-5 III 


T 




> 7.4 ± 0.5 


8 


4U 1543 - 475 


l''l2 


A2 V 


RT 




8.5 10.4 




XTE J1550 - 564 


1''55 


K3 III 


RT 


fiQ 


10.5 ± 1.0 




XTE J1650 - 500 


7''63 


K4 V 


T 




4.0 H- 7.3 


9 


GRO J1655 - 40 


2''62 


F3-6 IV 


RT 


mQ 


6.3 ± 0.5 




H 1705 - 250 


12'' 54 


K5 V 


T 




5.7 ^ 7.9 




GRO J1719 - 24 




MO-5 V 


T 




> 4.9 


10 


XTE J1859 + 226 


9'' 16 


- G 5 


T 




8 -MO 




GRS 1915 + 105 


33''5 


K-M III 


RT 




14 ± 4 




GS 2000 + 251 


8'' 3 


K5 V 


T 




7.1 ^ 7.8 




GS 2023 + 338 


6^*46 


KO IV 


RT 




10.0 ^ 13.4 





NOTES: 
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Port — orbital period 
Opt. Sp — optical spectrum 
X-R — X-ray variability 
C — comments 

Mbh— mass of black hole component 



Ref — references 

T — transient 

RT — recurrent transient 

pers — persistent 

yuQ — microquasar 



The errors or ranges for Mbh are in most cases quoted after original references. The detailed 
discussion of these estimates is given in Ziolkowski (2003). 
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(2005), (2) Hillwig and Gies (2006), (3) Casares et al. (2006), (4) Hynes et al. (2003), (5) Shahbaz 
et al. (2004), (6) Jonker et al. (2005), (7) Gelino et al. (2006), (8) Casares et al. (2004), (9) Orosz 
et al. (2004), (10) Masetti et al. (1996). 



2.2 Black Hole Candidates from Microlensing Events 

Among several hundreds of microlensing events observed so far there are about 30 so 
called "paralax events". These events are long enough to show the magnification fluctu- 
ations, reflecting the orbital motion of the Earth around the Sun. This effect permits to 
calculate the "microlensing parallax" which is a measure of the relative transverse mo- 
tion of the lens with respect to the observer. Assuming standard model of the Galactic 
velocity distribution, one is then able to perform a likelihood analysis, which permits 
to estimate the distance and the mass of the lens. With the help of the above analysis, 
three long events were selected as, possibly, caused by black hole lenses. The candi- 
dates are: MACHO-98-BLG-6 (probable mass of the lens ~ 3 13Mq, Bennett et al., 
2002a), MACHO-96-BLG-5 (probable mass of the lens ~ 3-=-16M0, Bennett et al., 2002a) 
and MACHO-99-BLG-22 = OGLE-1999-BUL-32 (probable mass of the lens ~ lOOM©, 
Bennett et al., 2002b). Only the last of them seems to be a robust candidate. I will also 
add, that Paczynski (2003) promised more BH lenses from OGLE III project in some 2 
-j- 3 years. OGLE III detects currently more than 500 events per year and, among them, 
some 20 30 paralax events. Based on the present (rather poor) statistics, we might 
expect that few of them (per year) should be BHCs. However, as no new firm detections 
were reported so far, it seems, that Paczynski's prediction was slightly overoptimistic. 

3 INTERMEDIATE MASS BLACK HOLES 

The case for the existence of intermediate mass black holes (IMBHs) is still not very 
strong but slowly it gets stronger. There are two classes of objects where we expect 
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IMBHs to be present: some ultraluminous compact X-ray sources (ULXs) and centers of 
some globular clusters (GCs). 



3.1 Ultraluminous X-Ray Sources (ULXs) 

The term ULXs is probably a sort of an umbrella covering several different classes of 
objects. However, the evidence is growing that one of these classes is, most likely, a class 
of XRBs containing IMBHs. 

Below, I briefly list the arguments supporting this hypothesis: 

• Some of them are confirmed XRBs 

• Some of these XRBs must be massive XRBs since they contain young massive 
optical components (08 V in M 81 X-1 (Liu et al., 2002), mid-B sg in M 101 X-1 (Kuntz, 
2005), BO lb in NGC 5204 X-1 (Liu et al., 2004)) 

• Many of the ULXs are found in star forming regions and young stellar clusters. 

• Many of them exhibits X-ray variability on time scales hours to years. 

• X-ray emission of some ULXs shows QPOs on time scales of few to few tens seconds, 
consistent with accretion discs around IMBHs. 

• X-ray spectra of many ULXs are consistent with relatively cool accretion discs 
around IMBHs. 

• Energy input into giant ionization nebulae surrounding many ULXs exclude sub- 
stantial beaming (Ho II X-1, M 81 X-9, M 81 X-6, NGC 1313 X-2, NGC 1313 X-1, NGC 
5408 X-1, IC 342 X-1, NGC 5204 X-1 and others). For a discussion of these objects see 
Pakull and Mirioni (2003) and Miller et al. (2005). 

3.1.1 Mass estimates for the compact components of ULX XRBs 

The mass of such objects can be estimated from the X-ray luminosities and the QPO 
frequencies (if observed). No estimate based on the mass function is available so far. 

3.1.2 M82 X-1 - the strongest ULX candidate for containing an IMBH 

The observed X-ray emission of this source corresponds to an isotropic luminosity of 
(2.4 16) X 10^" ergs/sec (which corresponds to the mass of 150^ lOOOM©, if the source 
emits at the Eddington level). The source exhibits 0.054 Hz and 0.114 Hz QPOs. The 
source is, most likely, an X-ray binary with an orbital period of ~ 62 days. It belongs to 
the young stellar cluster MGG-11 (7 to 12 Myr old). 

The analysis of these properties leads to the conclusion that, most likely, the system 
contains an IMBH of 200 to 5000 Mq accreting from a giant of ~ 25M0, filling its Roche 
lobe (Patnmo et al., 2006). 

The mass estimates for two other probable ULX IMBHs are given in Table 2. 
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3.2 Globular Clusters 

Modeling of the gravitational field in the central regions of some GCs indicates that they 
contain fairly massive 10^ to 10^ Mq) compact objects. It is likely that these objects 
are IMBHs (although, in some cases, a very dense cluster of neutron stars cannot be 
ruled out). 

3.2.1 Mass estimates for the compact objects in the centers of some GCs 

The principal method is the analysis of the brightness profiles of the central regions 
of GCs. Additional information might be obtained from the estimate of the velocity 
dispersion in the cores of the clusters. 

The useful parameter for detecting a probable presence of a central black hole during 
the preliminary analysis of the brightness profiles is the ratio of core radius to the half 
mass radius rc/vh- Trenti (2007) analyzed the dynamical evolution of a GC under a 
variety of initial conditions. She found, that for a cluster consisting initially from single 
stars only, the final (after relaxation) value of rc/rh was ~ 0.01, for a cluster containing 
10 % of binaries this value was ~ 0.1, but fot the cluster containing an IMBH the value 
of Tc/vh was ~ 0.3. These results confirmed earlier conclusions that IMBH clusters have 
expanded cores. Trenti considered subsequently 57 dynamically old (relaxed) GCs and 
found that for at least half of them the value of rc/rh is ^ 0.2, which implies the presence 
of an IMBH. It seems, therefore, that a substantial fraction of old GCs contains IMBHs. 
The case, however, is not closed, since a very dense cluster of neutron stars can mimic 
the gravitational potential of an IMBH. 

More detailed analysis of the brightness profiles leads to some quantitative estimates 
of the masses of the probable central BHs. Some of these estimates are given in Table 2. 

As far, as the velocity dispersion is concerned, it was noted (Gebhardt et al., 
2002), that GCs obey the relation (or, rather, an extension of it) between the velocity 
dispersion in the core and the mass of the central BH, found earlier for the galaxies. The 
relation might be useful for preliminary estimates of the central black hole mass, if the 
data about velocity dispersion are available. 

Finally, one should mention that other techniques were also used to estimate the 
central BHs masses in the stellar clusters. This includes the kinematics of milisecond 
pulsars (NGC 6752, Ferraro et al., 2003) and kinematics of the massive hot stars in the 
central region (IRS 13, Maillard et al., 2004). The results are also given in Table 2. 

3.3 An ULX in a Globular Cluster? 

At the end of this section, one should mention an object that might be an ULX in the 
GC. This object is a bright X-ray source in the unnamed GC which belongs to the Virgo 
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Cluster giant elliptical galaxy NGC 4472 (Maccarone et al. , 2007) . The source luminosity 
is ~ 4 X 10"^^ ergs/sec (which corresponds to the mass of ~ 25-^3OM0, if the source emits 
at the Eddington level). The source exhibits X-ray luminosity variability by a factor of 
7 in a few hours, which excludes the possibility that the object is several neutron stars 
superposed. It seems likely that the GC in question contains an ULX, which harbors a 
fairly massive BH (although, perhaps, not an IMBH yet). 

Tab. 2 — Masses of some intermediate mass black holes 



Name 


Mbh/Mq 


Ref 


Ultraluminous X-ray sources 


M 82 X-1 


200 ^ 5000 


1 


MCG 03-34-63 X-1 


^ 2000 


2 


Cartwheel N.IO 


^ 1000 


3 


Stellar clusters 


Gl 


- 20 000 


4 


M15 


~ 2000 


5 


u> Can 


~ 50 000 


6 


NGC 6752 


fewx (100 ^ 1000) 


7 


IRS 13 


~ 1300 


8 



REFERENCES: 
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et al. (2004). 



4 SUPERMASSIVE BLACK HOLES 



The term supermassivc black holes (SBHs) is used for BHs with the masses in the range 
~ 10^ to IQ^^Mq. Initially, such objects were believed to reside only in the centers of 



Masses of Black Holes in the Universe 9 

Active Galactic Nuclei (AGNs). During the last decade, the evidence was accumulated, 
which indicates, that virtually every "normal" galaxy contains in its center a SBH. 

4.1 Methods of Mass Estimates for Black Holes in the Centers of 
Galaxies 

At present, the four methods, listed below, are being used: 

• Kepler's law 

• Mbh — Mbuige relation 

• "Reverberation" method 

• "Variance" method 

4.2 Kepler's law 

This method is based on measuring the motions of the objects following the Keplerian 
orbits around the central BH. The objects might be either individual stars or stellar 
aggregates or water masers in the accretion disc. This method produces the most precise 
estimates, especially, if applied to individual stars or water masers. 

In particular, the mass of Sgr A*, the SBH in the center of our Galaxy was estimated 
from the motions of the several stars. In this case, both radial velocities and astrometric 
positions could be measured. These observations demonstrated that the common focus 
of the elliptical stellar orbits must harbor an invisible (except the rare X-ray or IR flares) 
mass of (3.6 ± 0.4) x 10^ ((Ghez et al., 2003). The pericenter distance of one of 
the stars (SO-16) is only about 8 light hours, so the size of the invisible object must be 
smaller. It could be nothing but BH. 

For our twin galaxy M31, the orbits of individual stars could not be used and only 
the kinematics of stellar aggregates was investigated. Prom the complex modeling of the 
gravitational potential, the value 1.4x10^ Mq was derived for the mass of the central 
BH. 

For a handful of nearby galaxies water maser sources were observed within the accre- 
tion disc (sometimes a few different sources at different distances from the center for a 
given galaxy). Their orbital motions yield quite precise estimate of the central mass. As 
an example, one may quote the galaxy NGC 4258 for which the value (3.9 ± 0.1) x 10^ 
Mq was found for the mass of the central BH (Herrnstein et al., 1999). 

4.3 Mbh — Mbuige relation 

Magorrian ct al. (1998) foimd an empirical relation between the mass of the bulge of the 
galaxy and the mass of the central BH. The more recent version of this relation may 
be found e.g. at Haring and Rix (2004). As the mass of the bulge is estimated from 
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the stellar velocities dispersion, this relation is, in fact, the relation between the stellar 
velocities dispersion and the mass of the central BH (which is obeyed also by some GCs, 
as was mentioned in the section 3.2.1). This relation may be quantified as: 

Mbh = 1-35 X 10^(a/200 km/sf -^^ Mq (7) 

4.4 "Reverberation" method 

This method can be applied in the case of active galaxies, where the previous method 
cannot be used (the emission from the active nucleus overshines the stellar component and 
the velocity dispersion which is based on stellar absorption lines cannot be determined). 
The reverberation method is essentially based on Kepler's law (Mbh = v'^R/G) applied, 
in this case, to broad line region (BLR). The velocity v of the matter in this region is 
estimated from the width of the H/3 emission line and the distance from the center R 
from the time delay of the variability of the emission from BLR with respect to that from 
the central region. Since both the radius R and the velocity v (or the H/3 width AA) were 
found empirically to scale with the total luminosity of the nucleus, the relevant formula 
can be simplified to (Kaspi et al., 2000): 

Mbh = 5.71 x 10^^44°-^^^ Mq (8) 

where L44 is the luminosity in the units of 10^^ erg/s. 

This method was successfully applied to many active galaxies. Unfortunately, its 
accuracy is rather low, since the result depends on the inclination of the plane of BLR 
(which is unknown) and on the assumption about the geometry of motions of the gas in 
BLR (which probably is not exactly keplerian). 

4.5 "Variance" method 

This method (Papadakis ct al., 2004; Nikolajuk et al., 2004) is based on the analysis of the 
variability in the X-ray band (long sequences of good quality observations arc necessary). 
The mass of the central object can be estimated from the high frequency break in the 
power density spectrum (the frequency scales approximately linearly with the mass. This 
method is completely independent of the previous one. The new calibration of the relevant 
formula (Nikolajuk et al., 2006) was based on the variability of the system Cyg X-1 (and 
on the new determination of its mass (Ziolkowski, 2005)). It is, certainly, encouraging 
that the masses determined by two completely independent methods (reverberation and 
variance) agree well with each other (Nikolajuk et al., 2006). 

Mass estimates were carried out for a large number of SBHs. Short selection of some 
results is given in Table 3. It is worth noticing that the range of masses covers more than 
five orders of magnitude. 
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Name 


Mbh/Mq 


Ref 


Low mass supermassive black holes 


NGC 4395 


3.6 X 10^ 


1 


NGC 4051 


~ 5 X 10^ 


2 


Sgr A* 


3.6 X 10® 


3 


High mass supermassive black holes 


[BH89] 1346-036 


0.9 X 10^° 


4 


LBQS 0109+0213 


1.0 X 10^° 


4 


[BH89] 0329-385 


1.3 X 10^° 


4 


2QZ J222006. 7-280324 


1.4 X 10^° 


4 


[BH89] 1246-057 


1.7 X 10^° 


4 


TON 618 


6.6 X 10^° 


4 



REFERENCES: 

(1) Peterson et al. (2005), (2) Shemmer et al. (2003), (3) Ghez et al., (2003), (4) Shemmer et 
al. (2004). 
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